Pseudo-diffusive magnetotransport in graphene 
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Transport properties through wide and short ballistic graphene junctions are studied in the pres- 
ence of arbitrary dopings and magnetic fields. No dependence on the magnetic field is observed at 
the Dirac point for any current cumulant, just as in a classical diffusive system, both in normal- 
graphene-normal and normal-graphene-superconductor junctions. This pseudo-diffusive regime is 
however extremely fragile respect to doping at finite fields. We identify the crossovers to a field- 
suppressed and a normal ballistic transport regime in the magnetic field - doping parameter space, 
and provide a physical interpretation of the phase diagram. Remarkably, pseudo-diffusive transport 
is recovered away from the Dirac point in resonance with Landau levels at high magnetic fields. 



Low energy excitations in a monolayer of carbon atoms 
arranged in a honeycomb lattice, known as a graphene 
sheet, have the remarkable peculiarity of being governed 
by the 2D massless Dirac equation, which is responsi- 
ble for a variety of exotic transport properties as com- 
pared to ordinary metals. Particularly striking is that 
for clean, undoped graphene the density of states is zero, 
but not so the conductivity, which remains of the order of 
the quantum unit e^/h [l|, Another intriguing fact is 
that a wide and short strip of undoped graphene exhibits 
"pseudo-diffusive" transport properties in the absence of 
electron-electron interactions and impurity scattering ^ . 
By pseudo-diffusive it is meant that transport properties 
are indistinguishable from those of a classical diffusive 
system. These include the full transport statistics (in 
particular the Fano factor F = 1/3 and the conductance 
G oc W/ L 3], where W is the width and L is the length 
of the graphene strip), the critical current [1] and I-V 
characteristics ^] in Josephson structures, as well as the 
relation of the normal metal-superconductor conductance 
to the normal transmissions Q . The same behaviour can 
be expected in bilayer graphene In fact, all of the 
above similarities can be explained by noting that at the 
Dirac point (i.e., for undoped graphene) transport occurs 
entirely via evanescent modes with a transmission that is 
equal to the diffusive transport theory result (evaluated 
&t kpl = 1, I = mean free path [131) without quantum 
corrections [1, H, [1], 
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Here ky is the transverse momentum of the channel. 



In 



diffusive systems the above relation holds independently 
of an externally applied magnetic field in the limit of 
many channels (classical limit), for which any quantum 
weak localization correction is negligible ^] . 

The question we raise here is: Does the diffusive be- 
haviour of ballistic graphene persist in the presence of a 
magnetic field? We will show that for zero doping the 
equivalence is preserved for any magnetic field. Remark- 
ably for a ballistic system, the applied magnetic field does 



not affect the transport statistics (for any current cumu- 
lant) at the Dirac point. For graphene with disorder this 
has rece ntly also been demonstrated at the conductivity 
level [§, [ToT [Tlj . At sufficiently strong magnetic fields an 
exponentially small chemical potential is enough to enter 
a field-suppressed transport regime. However, at reso- 
nance with the Landau levels (LLs) the pseudo-diffusive 
behaviour is recovered for all current cumulants. For 
even higher dopings one observes a final crossover to the 
ballistic magnetotransport regime, since clean graphene 
then resembles a ballistic normal metal. 

The magnetic field introduces a fundamental quantum- 
mechanical length scale known as the magnetic length 
— {^/\gB\)^/'^ . In complete absence of scattering, 
localized LLs are well formed and ballistic transport is 
suppressed. The only contribution to transport in this 
regime comes from resonant tunneling exactly at Lan- 
dau energies. In usual metals, this happens when the 
cyclotron diameter 2rc — 2l^gkp is smaller than the rel- 
evant scattering length (set by system size, disorder or 
temperature) . If 2rc is of the order of or larger than the 
scattering length, then delocalized states contribute to 
transport leading to Shubnikov-de Haas oscillations and 
the Quantum-Hall effect. For wide and short ballistic 
strips the relevant scattering scale is the strip length L, 
and scattering on lateral boundaries can be ignored. In 
graphene, the lowest LL lies precisely at the Dirac point. 
Besides, at this point kp — Q and thus r,, — Q inde- 
pendently of the magnetic field. Therefore, in contrast 
to the normal metal strip and to the high doping limit, 
at the Dirac point no delocalized bulk transport should 
take place for any magnetic field and resonant tunneling 
should be field independent. 

To confirm this hand-waving picture we ana- 
lyze theoretically magnetotransport effects through 
normal-graphene-normal (N) and normal-graphene- 
superconductor (NS) wide ballistic junctions at arbitrary 
dopings and magnetic fields. From an experimental point 
of view, transport properties of lightly doped graphene 
in contact with superconductors is currently being in- 
vestigated [111- Moreover, the properties of graphene 
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in stron g m aKnetic fields are also a subject of great 
interest [13l Il4l . Il5| . i n p articular in relation to weak 
(anti)localization [li, [13 [H, [H, [H- Here we con- 
sider a clean graphene sheet of width W (assumed to be 
the largest lengthscale in the system) in the y-direction 
through which transport occurs in the a;-direction. For 
X < —L/2 it is covered by a normal contact, and for 
x > L/2 it is covered either by a superconducting con- 
tact or a normal one. The central region is lightly doped, 
leading to a finite Fermi energy /i measured relative to the 
Dirac point, which can be varied by an external gate volt- 
age. The contact regions are modeled as heavily doped 
graphene, with Fermi energy /ic conveniently fixed to 
infinity with respect to both /i and the superconduct- 
ing gap A. The boundary conditions in the y-direction 
are irrelevant 0] for large aspect ratios W/L ^ 1. We 
choose periodic boundaries for simplicity. A constant ex- 
ternal magnetic field B is applied perpendicular to the 
graphene sheet. We assume the electrodes to be magnet- 
ically shielded, e.g. by covering them with materials with 
high-magnetic permeability. In the Landau gauge we can 
write the vector potential as ^ = (0, Bx, 0) for |a;| < L/2, 
and constant in the contact regions. This gauge is con- 
venient since the motions in x- and y-direction are un- 
coupled and ky remains a good quantum number. We 
neglect Zeeman splitting, so that the electron spin only 
enters as a degeneracy factor of 2 in the following cal- 
culation. Finally, we note that edge currents generally 
give a negligible contribution to transport in the W ^ L 
limit. 

We will compute the N and the NS (Andreev) inverse 
longitudinal resistivity, p^^j: ~ G{L/W), expressed in 
terms of the conductances 13, 61 
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and the shot noise using corresponding expressions in 
terms of the transmission for normal conducting contacts 
Tk [2l| . Note that the above expressions for the NS case 
are valid only if T^^ is left-right symmetric, which is not 
in general true in the presence of a magnetic field. In 
our particular setup it does indeed turn out to be sym- 
metric. The transmission through the central region is 
obtained by imposing current conservation at the inter- 
faces, which translates into continuity of the wavefunc- 
tion. In the chosen gauge the scattering problem is effec- 
tively one-dimensional, the transverse mode profile e**^"^ 
being the same in all regions, so we will only discuss the 
x-dependence of the wavefunctions from now on. 

The contact region eigenstates at energy e = 

hvp^ fc^ -|- ky — fic + P with respect to the central strip 

Dirac point are given by 
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(3) 



The spinor lives in the space of the two triangular sub- 
lattices that conform the graphene hexagonal lattice. 



s = ±1 is the 'valley' quantum number (for the degen- 
erate K and K' points), and Zk = exp[iarg(fc2; -t- iky)], 
which tends to Zk ~ Sign(fca;) when /ic — * oo. 

The spinor -^^^.ix) = [Ku,s{x)Afu,s{^)V for the 
central region is determined by the ID Dirac equation 
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with a = {x+dx) / \/2, = {x—dx)/V2, x = x/ls+kylB, 
A — sign e and rig = {lB\e\/hvp)'^ /2. Canonical relations 
[a, a+] = 1 are satisfied. The above equation corresponds 
to the K valley (s = 1), while the s = — 1 equation is ob- 
tained by swapping d and d^. Since the central region 
is bounded, no integrability condition must be met and 
the eigenspectrum of ([4]) is continuous. The usual LL 
solutions, which correspond to integer rig, are thus com- 
plemented by a larger family of divergent wavefunctions, 
typically localized around the interfaces x — ±L/2, with 
arbitrary > [l^]. At the Dirac point {n^ — 0) the 
components of the spinor are uncoupled, and the two 
eigenstates for 5 = 1 are ^"^^^1(2;) = [0, exp(— J^/2)]^ 

and *^/i(a;) = [exp{+x^ /2),0]^ . The s = -1 solution 
has interchanged spinor components. At finite energy 
(rie > 0) the solutions to Eq. (j4|) become 
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They have been expressed in terms of the even and 
odd (in x) solutions h'^°{x) of the Klein-Gordon equa- 

'^"^'^^ - nh'n^°\i) 
malized so that a+ = \A 



tion a'^ahn"' {x) = nhn'''{x) [the square of nor- 
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ric function iF,{a,b,z) = 1 + f ^ + |[2±1)£! 
Sn — sign{sin[7r(n + l/2)/2]}, these are 
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Imposing continuity for each {ky^e} at the inter- 
faces results in the following s— independent transmission 
probability 
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FIG. 1: (Color online) Inverse longitudinal resistivity in units 
of 4eV7r/i (a) and bulk Fano factor (b) for the N (left) and NS 
(right) junctions as a function of the (absolute value of the) 
Fermi energy (in units of hvp / L) . Different curves correspond 
to different values oi L/Ib (where Ib = \/h/eB), which ranges 
from zero for red curve to four for the dark blue one in steps 
of one. The same for the ratio Gns''/Gn""' in the inset. 



are expressed in terms of the wavefunctions at the bound- 
aries = (^7^ + K^B^- This gives the gen- 
eral transmissions Tk^ for arbitrary doping and magnetic 
field. It reproduces previously known results at -B = 
and non-zero doping Q as well as Eq. ^ for = 0. 

In Fig. ma) we plot as a function of the Fermi 
energy for increasing values of the ratio L/Ib oc ^/B. 
We recover the results obtained without magnetic field, 
namely that ctn and cns [where a — Pxxi^ — 0)1 tend to 
the known quantum-limited minimal conductivity value 
ie'^/nh at zero doping, whereas for \iJ.\L/hvF ^ 1 the 
slope of the asymptotes tends to O.SStt and 0.257r for the 
NS and N junction respectively Q. Remarkably, as we 
increase -B, all the curves remain unchanged at the 
Dirac point /i = 0. The Dirac-point Fano factor in Fig. 
[Ijb) is also unaffected by magnetic fields, and takes the 
classical diffusive value (1/3 for the N and 2/3 for the NS 
junction). This happens for any current cumulant, since 
at /i = the transmission given in Eq. ([5]) reduces to Eq. 
(HI) independently of B. 

However, for /i > the resistivities and the Fano fac- 
tors do depend on the magnetic field. In particular, for 
2rc < L (and above a certain critical value of L/Ib) 
transport can take place only at resonance with the LLs 
{fiL/hvp = \/2nL/lB), while for other dopings p~ 



-1 



suppressed like e (^Z'^) /2 fQj- ^jjg ]\[ junction and as 
^~(l/Ib) £qj. fT^Q ]\jg Q-j^Q j'^Q width of the resonances 
at the LL energies vanishes for 2rc ^ L, as we con- 
sider no disorder Remarkably, at these reso- 
nances for large fields coincides with the one at the Dirac 
point Ae'^/nh, a theoretical value that is usually associ- 
ated strictly with zero doping and that is interestingly at 
odds with some experimental findings 1]. In fact it can 
be analytically demonstrated that not only the conduc- 
tance but the whole pseudo-diffusive transport statistics 
is recovered at the resonances for high magnetic fields. 
Under this perspective, the field-independent resistivity 
at /i = can be understood as due to resonant transport 
through the zero-th LL that remains pinned at the Dirac 
point. The field suppressed regime is apparent for small 
but finite doping in Fig. [TJa), where strongly de- 
creases with increasing value L/Ib- Correspondingly, 
the bulk Fano factor reaches the tunneling limit value (1 
for the N and 2 for the NS junction) as transport gets 
suppressed [see Fig. [IJb)], in which limit the noise of the 
edge currents not considered here could be visible. In- 
creasing the Fermi energy further one enters the regime 
2rc > L, where is composed of two parts. The first 
part is linear in pLL/hvp, in agreement with the scaling 
with L behaviour of a ballistic conductor subject to a 
magnetic field (L independent conductance) . In particu- 
lar, for sufficiently high dopings, all curves in Fig. [Ij^a) 
become parallel and tend to the same (average) slope 
as the zero-field conductivity. The second contribution 
to p~x is an oscillating part, which for 2rc > L can no 
longer be explained by the resonance with LLs, since in 
that regime the effect of the boundaries is dominating the 
level structure in the central region. In fact, for 2rc ^ L 
the oscillations become equally spaced and are explained 
rather by a Fabry-Perot type effect, connected to reso- 
nant tunneling through the structure. 

In the inset of Fig. [DJa) the ratio G^'^'^'/G^"''^ is plot- 
ted as a function of p for the same values of L/Ib as 
in the main panel. At the Dirac point the ratio goes to 
one. At /i > the suppressed magnetotransport mani- 
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til doping reaches the ballistic threshold and the ratio 
starts growing again, finally reaching its asymptotic value 
0.38/0.25 ^ 1.52. As explained in Ref. |[6i], this value is 
expected in normal ballistic systems with Fermi wave- 
length mismatch. Note again here that, for sufficiently 
suppressed G^""^, the edge contribution neglected 
here will dominate transport. 

All the previous behaviours can be condensed in a 
quantitative way in the phase diagram shown in Fig. [21 
It contains three regions corresponding to the three dif- 
ferent transport regimes, namely: pseudo-diffusive (red), 
field-suppressed (blue) and ballistic (green), in the L/Ib 
and pL/hvp = kpL parameter space. The correspond- 
ing crossover lines between regions are solid (dashed) for 
the N (NS) junction (note that the background colours 
correspond to the boundaries of the N case). The bound- 
aries for the pseudo-diffusive region have been calculated 
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FIG. 2: (Color online) Phase diagram representing the 
crossovers from localized, pseudo-difltusive and ballistic trans- 
port regimes in the field-doping parameter plane. Solid 
(dashed) lines represent the boundaries for a N (NS) junc- 
tion. LLs are labeled by n. 

assuming a maximum deviation of ±10% with respect 
to the Dirac point conductivity Ae'^/irh. At low fields 
the width of the pseudo-diffusive window that brackets 
the Dirac point is roughly field independent, whereas for 
L/Ib > 1.8 for N (1.35 for NS) the window closes down 
as exp[— (L/?b)^/4]. Physically this means that at these 
higher fields the quasi-diffusive transport regime is ex- 
tremely fragile respect to doping, and an exponentially 
fast crossover to the field-suppressed (localized) regime 
takes place. The boundaries of the latter (blue region) 
were set by a crossover criterion p~]: < 0.1(4e^/7r/i). Its 
spiked shape is due to the peaked contributions to the 
field-suppressed p^^ discussed in the analysis of Fig. (TJ 
which are produced by resonant tunneling through LLs. 
When the magnetic field is increased, the positions of 
these peaks shift to higher dopings, converging on ra- 
dial lines with slope 1/V2n, while their width decreases 
exponentially. Above a certain value oi L/Ib the pseudo- 
diffusive regime is recovered and the resonances are thus 



coloured in red. The third region (green) is character- 
ized by oc L at fixed field and doping (which would 
correspond to radial lines in the phase diagram), and is 
therefore a ballistic transport regime. As expected from 
the arguments in the introduction, the boundary of the 
field-suppressed region closely follows the ballistic thresh- 
old 2rc = L. Finally, intermediate regions (white) are 
characterized by strongly oscillating conductivities. 

In conclusion, by computing the general transmission 
probabilities through short and wide graphene junctions, 
we have found that the transport properties at the Dirac 
point exactly match those of a classical diffusive system 
even in the presence of a magnetic field, which actually 
does not affect transport at all at zero doping. This be- 
haviour, which is associated to the existence of a zero-th 
LL pinned at the Dirac point, is however found to be ex- 
ponentially fragile respect to doping for high fields. By 
analyzing inverse longitudinal resistivity and higher cur- 
rent cumulants we have identified and interpreted the 
three distinct regimes that appear at finite magnetic 
fields and dopings, corresponding to pseudo-diffusive, 
field-suppressed and ballistic transport, and computed 
the phase diagram for the N and NS junctions in the rele- 
vant field-doping parameter space. Transport resonances 
at the LL energies are found in the field suppressed 
regime, with and all higher bulk current cumulants 
saturating to the pseudo-diffusive Dirac point values at 
high fields. The width of these resonances decreases ex- 
ponentially with magnetic field, although broadening due 
to disorder in real samples is expected, thus facilitating 
experimental observation. The reappearance of pseudo- 
diffusive transport at finite doping could shed light on 
the I/tt discrepancy between experiments and theoreti- 
cal results for the conductivity at the Dirac point. 

We thank the support of G. Schon in promoting this 
collaboration. This work benefited from the financial 
support of the European Community under the Marie 
Curie Research Training Networks and the ESR pro- 
gram. F. G. acknowledges funding from MEG (Spain) 
through grant FIS2005-05478-C02-01 and the European 
Union Contract 12881 (NEST). 



[1] K. S. Novoselov et al, Science 306, 666 (2004). 
[2] K. S. Novoselov et al, Proc. Nat. Acad. Sc. 102, 10451 
(2005). 

[3] J. Tworzydlo et al, Phys. Rev. Lett. 96, 246802 (2006). 
[4] M. Titov and C. W. J. Beenakker, Phys. Rev. B 74, 

041401(R) (2006). 
[5] J. C. Cuevas and A. Levy Yeyati, Phys. Rev. B 74, 

180501(R) (2006). 
[6] A. R. Akhmerov and C. W. J. Beenakker, Phys. Rev. B 

75, 045426 (2007). 
[7] I. Snyman and C. W. J. Beenakker, Phys. Rev. B 75, 

045322 (2007). 
[8] C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 
[9] P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, 74, 



235443 (2006). 

[10] J. A. Verges et al (2006), cond-mat/0610201. 

[11] V. P. Gusynin and S. G. Sharapov, Physical Review Let- 
ters 95, 146801 (2005). 

[12] H. B. Heersche et al, Nature 446, 56 (2007). 

[13] K. S. Novoselov et al. Nature 438, 197 (2005). 

[14] Y. B. Zhang et al. Nature 438, 201 (2005). 

[15] N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 
Rev. B 73, 125411 (2006). 

[16] S. Morozov et al, Phys. Rev. Lett. 97, 016801 (2006). 

[17] A. F. Morpurgo and F. Guinea, Phys. Rev. Lett. 97, 
196804 (2006). 

[18] E. McCann et al, Phys. Rev. Lett. 97, 146805 (2006). 
[19] X. Wu et al (2006), cond-mat/0611339. 



5 



[20] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, 
076602 (2007). 

[21] Y. M. Blanter and M. Biittiker, Phys. Rep. 336, 2 (2000). 
[22] A. D. Martino, L. Dell'Anna, and R. Egger, Phys. Rev. 

Lett. 98, 066802 (2007). 
[23] A. R. Akhmerov and C. W. J. Beenakker (2006), cond- 



mat/0612698. 

[24] It is noteworthy that the diffusive transport theory actu- 
ally assumes weak disorder kpl 3> 1. Hence Eq. ([ij is an 
extrapolation ol the theory to a dirty metal limit. 



